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Sudden  Death  of  Entanglement 

Ting  Yu1* *  and  ].  H.  Eberly2* 

A  new  development  in  the  dynamical  behavior  of  elementary  quantum  systems  is  the  surprising 
discovery  that  correlation  between  two  quantum  units  of  information  called  qubits  can  be  degraded 
by  environmental  noise  in  a  way  not  seen  previously  in  studies  of  dissipation.  This  new  route  for 
dissipation  attacks  quantum  entanglement,  the  essential  resource  for  quantum  information  as  well 
as  the  central  feature  in  the  Einstein-Podolsky-Rosen  so-called  paradox  and  in  discussions  of 
the  fate  of  Schrbdinger's  cat.  The  effect  has  been  labeled  ESD,  which  stands  for  early-stage 
disentanglement  or,  more  frequently,  entanglement  sudden  death.  We  review  recent 
progress  in  studies  focused  on  this  phenomenon. 


Quantum  entanglement  is  a  special  type  of 
correlation  that  can  be  shared  only  among 
quantum  systems.  It  has  been  the  focus 
of  foimdational  discussions  of  quantum 
mechanics  since  the  time  of  Schrodinger  (who 
gave  it  its  name)  and  the  famous  EPR  paper  of 
Einstein,  Podolsky,  and  Rosen  ( 1 ,  2).  The  degree 
of  correlation  available  with  entanglement  is  pre¬ 
dicted  to  be  stronger  as  well  as  qualitatively  dif¬ 
ferent  compared  with  that  of  any  other  known 
type  of  correlation.  Entanglement  may  also  be 
highly  nonlocal — e.g.,  shared  among  pairs  of 
atoms,  photons,  electrons,  etc.,  even  though  they 
may  be  remotely  located  and  not  interacting  with 
each  other.  These  features  have  recently  promoted 
the  study  of  entanglement  as  a  resource  that  we 
believe  will  eventually  find  use  in  new  approaches 
to  both  computation  and  communication,  for  ex¬ 
ample,  by  improving  previous  limits  on  speed 
and  security,  in  some  cases  dramatically  (3,  4). 

Quantum  and  classical  correlations  alike  al¬ 
ways  decay  as  a  result  of  noisy  backgrounds  and 
decollating  agents  that  reside  in  ambient  envi¬ 
ronments  (5),  so  the  degradation  of  entanglement 
shared  by  two  or  more  parties  is  unavoidable 
(6-9).  The  background  agents  with  which  we  are 
concerned  have  extremely  short  (effectively  zero) 
internal  correlation  times  themselves,  and  their 
action  leads  to  the  familiar  law  mandating  that 
after  each  successive  half-life  of  decay,  there  is 
still  half  of  the  prior  quantity  remaining,  so  that  a 
diminishing  fraction  always  remains. 

However,  a  theoretical  treatment  of  two-atom 
spontaneous  emission  (10)  shows  that  quantum 
entanglement  does  not  always  obey  the  half-life 
law.  Earlier  studies  of  two-party  entanglement  in 
different  model  fonns  also  pointed  to  this  fact 
(11-15).  The  tenn  now  used,  entanglement  sud¬ 
den  death  (ESD,  also  called  early-stage  disentan¬ 
glement),  refers  to  the  fact  that  even  a  very  weakly 
dissipative  environment  can  degrade  the  specif¬ 
ically  quantum  portion  of  the  correlation  to  zero 

department  of  Physics  and  Engineering  Physics,  Stevens  In¬ 
stitute  of  Technology,  Hoboken,  NJ  07030-5991,  USA.  Roch¬ 
ester  Theory  Center  and  Department  of  Physics  and  Astronomy, 
University  of  Rochester,  Rochester,  NY  14627-0171,  USA. 

*E-mail:  ting.yu@stevens.edu  (T.Y.),  eberly@pas.rochester.edu 
(J.H.E.) 


in  a  finite  time  (Fig.  1 ),  rather  than  by  successive 
halves.  We  will  use  the  tenn  “decoherence’ ’  to 
refer  to  the  loss  of  quantum  correlation,  i.e.,  loss 
of  entanglement. 

This  finite-time  dissipation  is  a  new  fonn  of 
decay  (16),  predicted  to  attack  only  quantum  en¬ 
tanglement,  and  not  previously  encountered  in  the 
dissipation  of  other  physical  correlations.  It  has 
been  foimd  in  numerous  theoretical  examinations 
to  occur  in  a  wide  variety  of  entanglements  in¬ 
volving  pairs  of  atomic,  photonic,  and  spin  qubits, 
continuous  Gaussian  states,  and  subsets  of  mul¬ 
tiple  qubits  and  spin  chains  (1 7).  ESD  has  already 
been  detected  in  the  laboratory  in  two  different 
contexts  (18, 19),  confirming  its  experimental  re¬ 
ality  and  supporting  its  universal  relevance  (20). 
However,  there  is  still  no  deep  understanding  of 
sudden  death  dynamics,  and  so  far  there  is  no 
generic  preventive  measure. 


How  Does  Entanglement  Decay? 

An  example  of  an  ESD  event  is  provided  by 
the  weakly  dissipative  process  of  spontaneous 
emission,  if  the  dissipation  is  “shared”  by  two 
atoms  (Fig.  1).  To  describe  this  we  need  a  suit¬ 
able  notation. 

The  pair  of  states  for  each  atom, 
sometimes  labeled  (+)  and  (— )  or  (1) 
and  (0),  are  quantum  analogs  of 
“bits”  of  classical  infonnation,  and 
hence  such  atoms  (or  any  quantum 
systems  with  just  two  states)  are 
called  quantum  bits  or  “qubits.”  Unlike 
classical  bits,  the  states  of  the  atoms 
have  the  quantum  ability  to  exist  in 
both  states  at  the  same  tune.  This  is 
the  kind  of  superposition  used  by 
Schrodinger  when  he  introduced  his 
famous  cat,  neither  dead  nor  alive 
but  both,  in  which  case  the  state  of 
his  cat  is  conveniently  coded  by  the 
bracket  (+  <=>  — ),  to  indicate  equal 
simultaneous  presence  of  the  oppo¬ 
site  +  and  -  conditions. 

This  bracket  notation  can  be  ex¬ 
tended  to  show  entanglement.  Sup¬ 
pose  we  have  two  opposing  conditions 
for  two  cats,  one  large  and  one  small, 


and  either  waking  (W)  or  sleeping  (S).  Entangle¬ 
ment  of  idealized  cats  could  be  denoted  with  a 
bracket  such  as  [(Ws)  <=>  (SW)],  where  we  have 
chosen  large  and  small  letters  to  distinguish  a  big 
cat  from  a  little  cat.  The  bracket  would  signal  via 
the  term  (Ws)  that  the  big  cat  is  awake  and  the 
little  cat  is  sleeping,  but  the  other  tenn  (Sw) 
signals  that  the  opposite  is  also  tme,  that  the  big 
cat  is  sleeping  and  the  little  cat  is  awake. 

One  can  see  die  essence  of  entanglement  here: 
If  we  leam  that  the  big  cat  is  awake,  the  (5w)  tenn 
must  be  discarded  as  incompatible  with  what  we 
learned  previously,  and  so  the  two-cat  state  re¬ 
duces  to  (Ws).  We  immediately  conclude  that  the 
little  cat  is  sleeping.  Thus,  knowledge  of  the  state 
of  one  of  the  cats  conveys  infonnation  about  the 
other  (21).  The  brackets  are  symbols  of  infor¬ 
mation  about  the  cats’  states,  and  do  not  belong  to 
one  cat  or  the  other.  The  brackets  belong  to  the 
reader,  who  can  make  predictions  based  on  the 
information  the  brackets  convey.  The  same  is  tme 
of  all  quantum  mechanical  wave  functions. 

Entanglement  can  be  more  complicated,  even 
for  idealized  cats.  In  such  cases,  a  two-party  joint 
state  must  be  represented  not  by  a  bracket  as 
above,  but  by  a  matrix,  called  a  density  matrix 
and  denoted  p  in  quantum  mechanics  [see  (22) 
and  Eq.  S3],  When  exposed  to  environmental 
noise,  the  density  matrix  p  will  change  in  time, 
becoming  degraded,  and  the  accompanying  change 
in  entanglement  can  be  tracked  with  a  quantum 
mechanical  variable  called  concurrence  (23),  which 
is  written  for  qubits  such  as  the  atoms  A  and  B  in 
Fig.  1  as 


C(p)  =  max[0,g(f)] 


(1) 


where  Q(t)  is  an  auxiliary  variable  defined  in 
tenns  of  entanglement  of  formation,  as  given  ex¬ 
plicitly  in  Eq.  S4.  C  =  0  means  no  entanglement 
and  is  achieved  whenever  Q(t)  <  0,  while  for 


Cavity  a 


Cavity  b 


Fig.  1.  Curves  show  ESD  as  one  of  two  routes  for  relaxation  of 
the  entanglement,  via  concurrence  C(p),  of  qubits  A  and  B  that 
are  located  in  separate  overdamped  cavities. 
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Fig.  2.  Atom-atom  entanglement  is  plotted  as  a  function  of  time  for 
k  values  in  the  range  0  to  1.  [Adapted  from  (20)]  For  all  values  of  k 
less  than  V3,  the  half-life  rule  is  obeyed,  but  for  k  between  Vi  and  1, 
it  is  not.  For  those  values,  the  curves  show  ESD,  i.e.,  becoming  zero 
in  a  finite  time  and  remaining  zero  thereafter.  The  two  curves 
marked  with  arrows  are  similar  to  the  curves  in  Fig.  1.  Time  is 
represented  byp  =  1  -  exp  (-rt). 


maximally  entangled  states  C  =  1,  and  C  is 
limited  to  the  range  1  >  C  >  0. 

In  the  case  of  spontaneous  emission  there  is 
no  environment  at  all,  except  for  the  vacuum.  The 
vacuum  can  still  have  a  noisy  degrading  effect 
through  its  quantum  fluctuations,  which  cannot 
be  avoided,  so  both  atoms  in  Fig.  1  must  eventually 
lose  their  excitation  and  come  to  their  ground  states 
with  the  rate  T.  Then  their  state  is  simply  ( — ),  a 
completely  disentangled  situation  because  learning 
that  die  state  of  one  atom  is  (— )  does  not  change 
our  information  about  the  other,  also  (-).  Thus, 
disentanglement  is  the  eventual  fate  of  the  pair. 

The  question  is,  how  quickly  do  they  meet 
their  fate?  For  the  initial  density  matrix  shown  in 
Eq.  S5,  the  answer  is  supplied  by  the  surface 
graphed  in  Fig.  2,  which  shows  possible  path¬ 
ways  for  entanglement  dissipation  as  a  function 
of  time.  The  k  axis  shows  that  the  time  evolution 
of  entanglement  depends  on  the  value  of  the  pa¬ 
rameter  that  encodes  the  initial  probability  for  the 
two  atoms  to  be  in  the  doubly  excited  state  (++). 
The  two  extreme  concurrence  curves  for  k  =  1 
and  k  =  0  are  the  ones  already  shown  in  Fig.  1. 

The  sudden  death  behavior  shown  in  the  right 
highlighted  curve  of  Fig.  2  is  a  new  feature  for 
physical  dissipation  (16,  20)  and  is  induced  by 
classical  as  well  as  quantum  noises  (24).  It  is 
counterintuitive,  based  on  all  previous  single-atom 
experience,  because  spontaneous  emission  is  a 
process  diat  obeys  the  half-life  mle  rigorously  for 
individual  atoms.  But  it  turns  out  that  the  two- 
qubit  correlation  does  not  follow  the  one-qubit 
pattern.  That  is,  the  sudden  death  does  not  come 
from  a  shorter  half-life;  the  entangled  joint  cor¬ 
relation  does  not  even  have  the  half-life  property. 

As  reported  in  (18),  the  first  experimental 
confirmation  of  ESD  was  made  with  an  all-optical 
approach  focusing  on  photonic  polarization.  It  was 
achieved  by  the  tomographic  reconstruction  of 
p (t),  and  from  it  the  Q(t)  variable,  and  thus  the 


concurrence  C(p).  In  the  exper¬ 
iment,  both  amplitude  and  phase 
noises  that  can  degrade  entangle¬ 
ment  were  realized  by  combin¬ 
ing  beam  splitters  and  mirrors. 

Can  Sudden  Death  Be  Avoided 
or  Delayed? 

The  issue  of  how  to  avoid  ESD- 
type  decorrelation  in  a  realistic 
physical  system  is  incompletely 
resolved  at  this  time.  A  number 
of  methods  are  known  to  pro¬ 
vide  protection  against  previous¬ 
ly  known  types  of  decorrelation 
(3).  Some  methods  have  clas¬ 
sical  analogs  in  information  the¬ 
ory.  One  engages  appropriately 
designed  redundancy  and  is  known 
as  quantum  error  correction.  An¬ 
other  relies  on  using  a  symmetry 
that  can  isolate  entanglement 
from  noise,  effectively  providing 
a  decoherence-free  subspace  to 
manage  qubit  evolution. 

Error  correction  is  most  useful  when  the  dis¬ 
turbing  noise  is  below  some  threshold  (25).  In 
practice  error  correction  can  be  complicated,  be¬ 
cause  a  noisy  channel  is  a  dynamical  process  and 
its  physical  features  are  often  not  fully  under¬ 
stood  or  predictable.  An  example  is  atmospheric 
turbulence  during  open-air  communication.  An¬ 
other  issue  is  the  cost  associated  with  providing 
redundancy.  Additionally,  it  has  been  reported  that 
some  quantum  error  correction  algorithms  could 
actually  promote  rather  than  mitigate  ESD  (26). 
Symmetries  that  avoid  decoherence  by  providing 
isolation  from  noise  during  evolution  (27)  have 
also  been  examined  as  a  way  to  postpone  or  avoid 
ESD,  but  knowledge  of  the  noise  to  be  combatted 
appears  unlikely  to  be  available  because  qubits 
remote  from  each  other  would  rarely  share  either 
symmetry  properties  or  noise  descriptions. 

Other  methods  considered  use  dynamic  ma¬ 
nipulation  such  as  mode  modulation  (28)  or  the 
quantum  Zeno  effect  (29),  which  can  be  regarded 
as  extensions  of  the  so-called  bang-bang  method 

(30) .  Another  proposal  is  to  use  feedback  control 

(31)  to  prevent  ESD  and  other  decoherence 
effects  in  the  presence  of  hostile  noise.  None  of 
these  methods  is  perfect,  but  they  can  be  more 
effective  if  designed  for  specific  noise  avoidance. 

Does  the  Number  of  Noises  Matter? 

Nonlocal  entanglement  raises  the  issue  of  ESD 
triggered  by  different  noise  processes  and  can  refer 
to  more  than  one  noise  source  acting  together  on 
colocated  entangled  qubits,  or  to  independent  noise 
sources  acting  separately  on  remotely  located  mem¬ 
bers  of  a  qubit  pair.  The  rate  of  dissipation  in  the 
presence  of  several  noise  sources  is  normally  the 
sum  of  the  individual  dissipation  rates.  More  ex¬ 
plicitly,  if  decay  rates  F  {  and  T2  come  from  the 
action  of  two  distinct  weak  noises,  then  when  the 
two  noises  are  applied  together  to  a  physical  sys- 
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tem,  the  resulting  relaxation  rate  is  simply  given 
by  the  sum  of  the  separate  rates:  F,  +  T2. 

However,  such  a  long-standing  result  does 
not  hold  for  entanglement  decay.  This  was  dis¬ 
covered  (16)  by  examining  entanglement  evo¬ 
lution  of  a  set  of  X-form  mixed-state  matrices 
(Eq.  S3)  with  d  =  0,  where  d  is  the  probability  that 
both  qubits  are  in  their  ground  states.  Straight¬ 
forward  calculations  for  the  entire  class  are  shown 
by  the  diagrams  in  Fig.  3,  illustrating  two  qubits 
exposed  together  to  amplitude  noise  as  well  as 
phase  noise.  The  top  two  time-dependent  curves 
show  that  the  application  of  either  noise  separately 
allows  long-running  entanglement  decay  of  the 
half-life  type  (no  ESD).  By  contrast,  the  bottom 
curve  reaches  zero  in  a  finite  time.  That  is,  the  com¬ 
bined  effect  of  the  two  noises  causes  ESD.  The 
caption  explains  the  colored  squares.  This  illustrates 
the  “supervulnerability"  of  paired-qubit  entangle¬ 
ment  when  attacked  by  different  independent 
noises.  This  result  is  universal  in  the  sense  that  it 
continues  to  hold  (16)  even  if  the  two  noises 
attack  one  of  the  qubits  but  not  the  other,  and  also 
if  each  of  the  two  qubits  is  remotely  attacked  by 
just  one  of  the  noises. 

Is  There  an  "Anti-ESD"  or  Rebirth  Effect? 

Special  circumstances  are  needed  to  see  “anti- 
ESD,"  the  creation  or  rebirth  of  entanglement 
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Fig.  3.  The  top  two  time-dependent  curves  show 
exponential  (smooth  half-life  type)  decay  of  concur¬ 
rence  for  a  qubit  pair  exposed  to  phase  noise  and 
to  amplitude  noise,  respectively.  The  bottom  curve 
shows  nonsmooth  decay,  i.e.,  ESD  occurs  for  the 
qubit  pair  when  both  noises  are  acting  together. 
The  color-coded  squares  apply  to  any  two-party  X 
matrix  (Eq.  S3)  having  d  =  0.  They  present  the 
predicted  results  for  the  entire  accessible  physical 
domain,  which  means  throughout  1  >  a  >  0  and 
1  >  Izl  >  0.  The  blue  zones  labeled  EXP  designate 
domains  where  smooth  exponential  (half-life  type) 
evolution  occurs,  and  the  orange  zones  show  where 
ESD  occurs.  The  left  and  middle  squares  apply  when 
amplitude  and  phase  noises  are  applied  separately, 
and  have  smooth  decay  regions,  while  the  orange 
square  at  right  shows  that  sudden  death  is  universal 
in  the  entire  region  for  anyX  matrix  (Eq.  S3). 
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from  disentangled  states.  By  imposing  the  right 
interactions  almost  anything  can  be  made  to 
happen,  but  we  are  concerned  with  evolution  of 
joint  information  in  a  pure  sense  and  focus  on 
two-party  entanglement  that  evolves  without 
mutual  interaction  or  communication. 

The  same  two-atom  situation  shown  in  Fig.  1 
can  be  made  relevant  to  anti-ESD.  In  solving  for 
the  surface  of  solutions  plotted  in  Fig.  2,  the 
cavities  were  taken  as  hilly  overdamped,  so  that 
any  photon  emitted  by  either  atom  was  irreversibly 
absorbed  by  the  walls,  but  they  could  also  be 
treated  as  undamped  mirror-like  cavities,  such  as 
used  in  the  Jaynes-Cummings  (JC)  model  for  light- 
atom  interactions  (32).  This  situation  produces  a 
periodic  sequence  of  perfect  rebirths  of  atom  AB 
entanglement  (33,  34).  An  early  mathematical 
model  of  two-qubit  evolution  (11)  can  be  inter¬ 
preted  as  treating  an  underdamped  cavity  and  also 
shows  rebirths. 

The  panels  in  the  top  row  of  Fig.  4  show 
rebirth  scenarios.  They  occur  for  states  that  are 
initially  of  the  cat  type,  such  that  both  atoms  are 
excited  and  both  are  un-excited  at  the  same  time. 
The  cat-type  bracket  for  them  is 

<t>a  =  [(++)cos  a  o  ( — )sin  a]  (2) 

where  different  values  of  sine  and  cosine  produce 
the  different  curves  in  the  figure  (35). 

Starting  from  the  photon  vacuum  state  in  each 
cavity,  the  JC-type  evolution  will  permit  only  zero 
or  one  photon  to  reside  in  each  cavity  at  any  later 
time,  so  each  of  the  two  modes  is  a  two-state 
system  (a  qubit),  and  counting  the  two  atoms  there 
are  now  four  qubits  on  hand.  This  provides  six 
concurrences  that  can  be  computed:  C48,  Cah, 
CAb,  CaB,  CAa,  and  CBh,  where  the  capital  letters 
identify  atoms  and  the  small  letters  identify  the 
photons  in  cavity  modes  a  or  b.  Concurrence  is 
defined  only  for  pairs  of  qubits,  not  quartets  of 
them,  so  the  label  CAB  implies  that  the  a  and  b 
degrees  of  freedom  are  not  available  to  observa¬ 
tion  and  have  been  ignored  (technically,  have  been 
traced  out),  and  for  photon  concurrence  Cab  the 
atomic  A  and  B  properties  are  traced  out,  and  so  on. 

This  idealized  model  provides  a  convenient 
framework  to  analyze  entanglement  in  a  simpli¬ 
fied  but  still  multi-qubit  framework.  As  shown  in 
the  top  half  of  Fig.  4,  ESD  takes  place  for  atom 
concurrence  (in  the  left  panel,  almost  all  C  -Ls  curves 
reach  zero  and  remain  zero  for  finite  times).  How¬ 
ever,  in  the  right  panel,  the  photon  concurrence 
Cab  behaves  in  a  manner  opposite  to  that  of  C'48, 
showing  anti-ESD.  That  is,  initially  Cah  is  zero, 
but  it  immediately  begins  to  grow.  The  photons 
jointly  experience  entanglement  “sudden  birth,” 
but  this  is  followed  by  ESD  a  half  cycle  later.  All  of 
this  occurs  via  pure  “informatics,”  i.e.,  without  en¬ 
ergy  exchange  or  other  interaction  between  the  sites. 

The  reason  for  the  rebirths  is  obvious — the 
photons  emitted  cannot  really  get  lost  among  the 
few  joint  states  available.  If  a  larger  number  of 
cavity  modes  would  be  provided,  a  longer  time 
would  be  needed  for  a  rebirth  to  be  complete,  and 


as  a  limiting  case,  the  cavities  producing  the  curves 
in  F  ig.  2  have  an  infinite  number  of  modes,  so  the 
lost  quantum  correlation  cannot  be  reborn  in  any 
finite  tune.  If  there  are  enough  states  available  in 
one  mode,  as  is  the  case  for  coherent-state  mode 
preparation,  then  ESD  and  true  long-time  reviv¬ 
als  are  also  predicted  (36). 

What  Are  the  Future  Prospects? 

Quantum  memory  banks.  Clearly,  ESD  can  be 
largely  ignored,  to  a  first  approximation,  when 
desirable  quantum  operations  can  be  manipulated 
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at  sufficiently  high  speed.  The  key  goal  of  mem¬ 
ory  is  opposite  to  that  of  speed,  i.e.,  to  preserve 
quantum  state  features  semi-indefinitely.  Quan¬ 
tum  memory  networks  will  be  sensitive  to  the 
consequences  if  ESD  occurs.  ESD  will  probably 
have  to  be  taken  into  account  if  practical  versions 
of  quantum  memories  are  built  to  operate  in 
mixed-state  configurations. 

Disentanglement  control.  Over  a  given  noisy 
channel,  it  appears  that  some  entangled  states 
may  be  more  robust  against  the  influence  of  noise 
than  others.  To  control  decoherence  optimally,  it 
will  be  useful  to  leant  how  to  identify  the  robust 
states  separately  from  the  fragile  ones.  Control 


issues  also  include  the  use  of  external  fields  to 
manipulate  qubit  states  (37)  as  in  gate  operations, 
to  create  transient  decoherence-free  subspaces,  as 
mentioned  already.  A  qualitatively  different  route 
to  combat  decoherence  specifically  of  the  ESD 
type  is  illustrated  in  Fig.  2.  Two  evolution  tracks 
are  highlighted  to  show  that  for  qubits  prepared 
with  the  same  value  of  initial  entanglement,  their 
concurrences  may  evolve  very  differently.  In  that 
illustration,  the  right  track  is  subject  to  ESD,  while 
the  left  one  is  not.  Decoherence  per  se  is  not 
avoided,  because  the  non-ESD  track  shows  steady 
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dissipation,  but  it  always  remains  finite.  Other 
examples  of  this  are  known,  and  in  some  cases  a 
purely  local  operation  (i.e.,  a  manipulation  of  only 
one  of  the  two  entangled  qubits)  can  be  under¬ 
taken  to  change  the  state  matrix  p  without  chang¬ 
ing  its  degree  of  entanglement,  but  in  a  way  that 
switches  the  evolution  trajectory  from  ESD  to 
non-ESD  (38),  effectively  putting  it  on  a  half-life 
track  as  in  the  figure.  Similar  studies  (39)  have 
examined  the  effect  of  local  operations  at  inter¬ 
mediate  stages  of  evolution.  Use  of  this  method 
requires  detailed  knowledge  of  the  state  matrix  p, 
which  may  not  be  practical,  particularly  at  times 
late  in  the  evolution. 


Fig.  4.  Entanglement  birth,  death,  and  rebirth.  [Adapted  from  (35)]  In  the  bottom  pair  of  panels  the  rise 
and  fall  of  AB  entanglement  is  exactly  compensated  by  the  fall  and  rise  of  ab  entanglement.  This  is  not  the 
case  in  the  top  pair  of  panels,  but  a  more  subtle  form  of  compensation  still  occurs,  as  reported  in  (35).  It 
involves  not  concurrence  but  the  auxiliary  variable  0(0  defined  in  Eqs.  1  and  S4. 
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Entanglement  invariants.  Entanglement  flow 
in  small  reservoirs  has  led  to  the  recent  dis¬ 
covery  of  entanglement  “invariants”  (55)  by  in¬ 
spection  of  the  curves  in  the  bottom  row  of 
Fig.  4,  which  repeats  the  top  row  except  that  a 
slightly  different  cat  state  is  used  for  the  atoms 
initially: 

'Fa  =  [(+- )cos  a  o  (-+)sin  a]  (3) 

In  the  bottom  left  panel  of  Fig.  4,  it  appears 
that  each  AB  atomic  concurrence  curve  is  com¬ 
pensated  at  all  times  by  the  corresponding  ab 
photonic  concurrence  curve  in  the  right  panel, 
one  going  up  as  the  other  falls. 

In  fact,  exact  compensation  can  be  confirmed 
analytically  (35),  but  one  notes  that  the  same  be¬ 
havior  does  not  appear  in  the  top-row  curves  in 
Fig.  4,  where  there  is  no  perfect  compensation  of 
ab  for  AB.  For  example,  it  is  easy  to  see  that  the 
AB  and  ab  red  curves  can  be  zero  at  the  same 
time.  A  natural  question  is,  where  does  the  miss¬ 
ing  information  go  in  that  case?  Because  the 
two-site  JC  model  is  unitary,  preservation  of  all 
four-qubit  information  is  guaranteed,  so  it  should 
be  located  “somewhere.”  Careful  examination 
shows  that  concurrence  is  not  conserved,  but 
rather  Q(t)  is  conserved,  spread  among  all  six  dif¬ 
ferent  types  in  that  case  (35). 

This  identification  of  four-particle  memory 
flow  channels  is  unusual  and  clearly  deserves 
future  examination  (40,  41).  One  can  say  about 
these  invariants  that  they  emerge  only  from  a 
kind  of  analytic  continuation  of  the  bipartite  con¬ 
currence  function  C(t)  to  un-physically  negative 
values,  which  is  pennitted  via  Q(t).  The  entan¬ 
glement  flow  issue  (42)  is  also  related  to,  and 
appears  to  expand  considerably,  the  theory  asso¬ 
ciated  with  entanglement  swapping,  which  is  under 
active  exploration,  and  has  been  realized  with 
particle  pairs  from  independent  sources  (43). 

Non-Markovian  noises.  Dissipative  entangle¬ 
ment  evolution  is  critically  dependent  on  the  types 
of  the  noises  acting  on  the  system.  Markov  envi¬ 
ronments  are  those  for  which  a  noise  signal  has 
no  self-correlation  over  any  time  interval,  and 
under  Markov  conditions  noise  typically  results 
in  a  quantum  irreversible  process.  Non-Markovian 
noise  arising  from  a  structured  environment  or 
from  strong  coupling  appears  to  be  more  fun¬ 
damental  (44, 45).  Recent  studies  have  suggested 
that  correlated  noises  may  cause  new  difficulties 
in  using  quantum  error  correction  codes  (46)  and 
dynamic  decoupling  (47).  Although  some  progress 
has  been  made,  extending  the  current  research  on 
ESD  into  physically  relevant  non-Markovian  sit¬ 
uations  remains  a  challenge.  ffigh-Q  cavity  QED 
and  quantum  dot  systems  are  two  possible  exper¬ 
imental  venues. 

Qutrits  and  beyond.  Many-qubit  entanglement 
and  entanglement  of  quantum  systems  that  are  not 
qubits  (i.e.,  those  having  more  than  two  states)  is 
largely  an  open  question,  and  one  that  is  embar¬ 
rassing  in  a  sense,  because  the  question  has  been 
open  since  quantum  mechanics  was  invented  in 
the  1920s.  There  is  still  no  known  finite  algo¬ 


rithm  for  answering  the  simple-seeming  question 
of  whether  a  given  mixed  state  is  entangled  or 
not,  if  it  refers  to  more  than  two  systems,  and  it  is 
answerable  for  two  mixed-state  systems  only  in 
the  case  of  pairs  of  qubits  (as  we  have  been  dis¬ 
cussing)  and  the  case  of  one  qubit  and  one  qutrit 
(a  three-state  system  such  as  spin-1).  Investiga¬ 
tion  into  ESD  of  qubit-qutrit  systems  has  begun 
(48,  49),  but  generalizations  to  many-qubit  sys¬ 
tems  are  daunting  tasks  due  to  both  technical  and 
conceptual  difficulties  (7-9). 

Topological  approach.  A- party  entanglement 
dynamics  will  presumably  become  simpler  to  pre¬ 
dict  if  a  computable  entanglement  measure  for  a 
mixed  state  of  more  than  two  qubits  can  be  dis¬ 
covered.  However,  an  alternative  approach  is  to 
avoid  dynamics  through  topological  analysis  (50). 
One  now  knows  (51,  52)  that  ESD  is  necessary 
(i.e.,  must  occur)  in  arbitrary  A-party  systems  of 
noninteracting  qubits  if  they  are  exposed  to  ther¬ 
mal  noise  at  any  finite  T  >  0  temperature.  The 
steady  state  of  any  noninteracting  A-qubit  system 
has  a  neighborhood  in  which  every  state  is  sepa¬ 
rable.  In  this  case,  any  prearranged  subsystem 
entanglement  will  inevitably  be  destroyed  in  a 
finite  time.  This  is  a  universal  result  showing  how 
entanglement  evolves  in  the  absence  of  external 
noise  control. 

Clearly,  the  holy  grail  for  research  on  entan¬ 
glement  dynamics  is  to  find  an  efficient  real-time 
technique  for  tracking  and  controlling  the  entan¬ 
glement  evolution  of  a  generic  many-qubit  system. 
Another  important  open  question  is  to  determine 
a  generic  method  for  direct  experimental  registra¬ 
tion  of  entanglement,  for  which  there  is  no  current 
answer.  We  believe  that  many  surprising  results 
are  awaiting  discovery. 
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We  report  calculations  of  the  long-time  evolution  of  quantum  entanglement.  We  employ  a  simple 
model  that  permits  partially  analytic  examination  of  entanglement  transfer  among  a  combination  of 
discrete-state  and  continuous-state  non-localized  degrees  of  freedom  of  a  four-part,  quantum  system. 
Among  the  results  are  the  development  of  a  novel  approximation  scheme,  the  re-appearance  of 
entanglement  after  lengthy  intervals  of  zero  entanglement,  and  formulas  for  relatively  slow  loss  of 
entanglement  and  duration  of  revival  intervals. 

PACS  numbers:  42.50.Ex,  42.50.Pq,  03.65.Ud,  42.5Q.dt. 


I.  INTRODUCTION 


Almost  75  years  after  its  description  by  Schrodinger  [1] 
there  is  still  no  general  understanding  of  entanglement 
within  quantum  theory  in  the  sense  that  measures  of  it 
are  still  being  developed.  Only  two-particle  entanglement 
is  well  quantified,  and  even  in  that  case  the  number  of 
states  that  can  be  taken  into  account  is  extremely  small 
unless  the  joint  state  is  pure.  The  behavior  of  few-state 
entangled  quantum  systems  takes  a  wide  variety  of  forms 
in  different  environments,  and  it  is  well  known  that  some 
of  them  are  being  studied  to  achieve  advances  in  com¬ 
puting  and  communications.  Whether  such  goals  are  of 
interest  or  not,  the  controlled  manipulation  of  entangle¬ 
ment  with  external  agents  is  one  approach  to  an  ultimate 
form  of  quantum  control,  and  it  remains  a  fascinating 
challenge. 

We  will  examine  aspects  of  what  we  can  call  pure  stor¬ 
age  of  quantum  information,  a  context  that  might  have 
relevance  in  computing  and  communication.  More  gen¬ 
erally  this  means  that  we  are  interested  in  evolution  of 
two-party  entanglement  of  simple  examples  of  quantum 
systems  such  as  two-level  atoms  and  coherent  states  of 
radiation,  but  in  this  case  under  the  conditions  of  mutual 
isolation  and  non-interaction  between  the  joint  informa¬ 
tion  storing  systems.  The  term  qubit  serves  as  a  use¬ 
ful  general  term  for  two-level  quantum  systems  whether 
in  atoms,  molecules,  impurity  sites  in  crystals,  quantum 
dots,  SQUIDs,  photon  polarizations  or  of  course  spins,  so 
we  will  adopt  it  for  a  discussion  of  long-time  t.wo-qubit 
evolution  under  the  no-interaction  stipulation. 

Earlier  discussions  usually  allowed  or  relied  on  mutual 
interactions  to  produce  entanglement  dynamics  or  con¬ 
trol  it.  Ivim,  et.  ah,  showed  [2]  that  an  incoherent  thermal 
field  can  create  entanglement  between  two  such  qubits. 
Entanglement  transfer  between  two  qubits  and  two  sep¬ 
arate  fields  was  examined  by  Zhou  and  Wang  [3],  In 
their  treatment  the  quantum  field  was  weak  rather  than 
strong.  The  evolution  of  entanglement  in  a  qubit-field 
system,  where  the  qubit  and  the  field  start  from  mixed 
states  was  examined  by  Ren  dell  and  Rajagopal  [4].  They 
aimed  to  calculate  the  entanglement  embedded  in  the  full 
system,  and  because  of  the  lack  of  an  entanglement  mea¬ 


sure  for  2xoo  systems  they  calculated  a  lower  bound  for 
the  concurrence  instead.  The  controlled  manipulation  of 
entanglement  with  external  agents  is  a  common  goal  of 
most  of  these  studies. 

We  first  review  entanglement  dynamics  in  an  ideal¬ 
ized  situation  in  order  to  establish  notation  for  control 
parameters,  time  scales,  etc.  Thus,  each  site  of  our  net¬ 
work  contains  a  qubit  in  the  physical  form  of  a  two-level 
atom  (ion  qubit,  spin  qubit  and  photonic  qubit  analogs 
are  easily  imagined) .  We  restrict  attention  to  the  dynam¬ 
ics  of  entanglement,  between  two  of  the  sites,  say  A  and 
B.  We  suppose  that  there  is  quantum  memory  present, 
in  the  qubits  that,  is  of  interest,  i.e.,  that  the  two  atoms 
have  been  entangled  in  some  way  before  being  inserted 
into  their  respective  cavities.  Field  modes  a  and  b  are 
available  at  the  A-B  sites  that  can  be  used  for  interact¬ 
ing  with  the  atoms  (i.e.,  for  “controlling”  them  and  their 
entanglement  externally).  This  scenario  is  suggested  in 
Fig.  1. 


Atom  A  Atom  B 


FIG.  1:  A  diagram  illustrating  our  system.  The  atoms  A  and 
B  are  placed  in  their  respective  cavities  a  and  b  after  being 
entangled.  Only  A-a  and  B-b  interactions  take  place  there¬ 
after,  causing  non-interacting  A-B  entanglement  to  change  in 
time. 
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II.  HAMILTONIAN  EQUATIONS 

The  Jaynes-Cummings  [5]  interaction  governs  evolu¬ 
tion  in  our  system.  The  Hamiltonian  (with  ft  =  1)  is 
given  by: 

Htot  =  y-rl4  +  g{a^<r^  +  a+a)  +  uia^a 

+Ycrf  +g{^(T-  +<r+b)  +ubH,  (l) 

where  u>o  is  the  transition  frequency  between  the  two  lev¬ 
els  of  the  atoms,  g  is  the  constant  of  coupling  between  the 
atoms  and  the  fields  and  ui  is  the  angular  frequency  of  the 
single-mode  field.  The  usual  Pauli  matrices  crf'B,  cr+’B 
and  describe  the  atoms,  while  at,  a  and  ftt,  b  are 

the  raising  and  lowering  operators  for  the  fields. 

Note  that  the  Hamiltonian  does  not,  support  any  in¬ 
teraction  between  atom  A  and  atom  B  or  between  mode 
a  and  mode  b.  Physical  realization  of  our  scenario  does 
not  appear  out  of  the  question,  as  the  Jaynes-Cummings 
model  has  been  realized  in  the  laboratory  in  several  well- 
known  ways  [6-9]. 

To  illustrate  the  approach  to  be  followed  when  the 
fields  are  taken  as  coherent  states,  we  start  with  a  sim¬ 
pler  example,  in  which  the  field  modes  are  initially  in 
their  vacuum  states:  |0a)  (I)  |0&),  and  the  two  atoms  are 
in  a  superposition  of  Bell  states.  For  the  vacuum  case 
the  coupling  terms  in  our  Hamiltonian  indicate  that  for 
creation  of  a  photon  in  a  cavity  the  atom  in  that  cav¬ 
ity  must  decay  to  a  lower  state.  Since  the  cavities  each 
contain  only  one  atom,  there  can  be  no  more  than  one 
photon  in  each  at  any  time.  This  means  that  the  cavities 
are  also  two-state  quantum  systems,  i.e. ,  qubits. 

The  eigenstates  of  the  JC  Hamiltonian  are  well  known. 
We  will  write  for  either  Aa  or  Bb 

H3c\4$)  =  Xt\4$).  (2) 


— it 
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FIG.  2:  Sketch  indicating  non-interacting  qubits  in  a  quantum 
storage  network.  Quantum  memory  in  the  form  of  two-qubit 
entanglement  between  any  two  relatively  remote  sites,  indi¬ 
cated  here  by  a  dashed  line,  can  be  altered  by  fields  local  to 
the  sites. 


We  will  denote  the  excited  and  ground  atomic  states  by 
|e)  and  | g),  and  denote  the  cavity  modes’  photon  states 
by  the  photon  number  n.  Then  the  JC  eigenvalues  are 


given  by: 

A±  =nw  +  i(A±  ^A2  +  Gl),  (3) 

where  A  =  u>o  —  w  is  the  detuning  and 

Gn  =  2g\/n,  (4) 

is  the  n-photon  Rabi  frequency. 

The  JC  eigenstates  have  the  following  form  as  super¬ 
positions  of  the  bare  atom  and  cavity  states: 

iV’a)  =  |3;0>  (5) 

H’i)  =  Cn|e;  n  —  1}  +  sn\g;  n)  (n  >  0)  (6) 

\4’n)  =  sn\e\n  —  1)  +  \c„\g;n)  (n  >  0).  (7) 


In  these  equations  we  have  introduced  some  convenient 
abbreviations: 


cn  =  cos($„/2)  and  %  =  sin(#„/2),  (8) 


where  the  rotation  angle  9n  can  be  identified  with  the 
Bloch  sphere  polar  angle  and  is  defined  in  the  usual  way: 


A 


COS  Un  = 


>/A2  +  G* 


and 


G„ 


•  (9) 


In  this  preliminary  example  we  will  need  the  true 
ground  state  |t/’o)  =  | g\  0)  and  the  two  dressed  states  for 
n  =  1.  These  three  states  are  closed  under  the  JC  Hamil¬ 
tonian  for  each  site.  In  other  words,  we  use  only  n  =  1 
in  the  equations  above,  so  the  subscript  n  can  mostly  be 
ignored  and  we  will  frequently  drop  it  (Xn  — >  A,  cn  — >  c, 
etc.). 


III.  MEASURE  OF  PAIRWISE 
ENTANGLEMENT 

The  JC  dressed  states  are  atom-photon  locally- 
entangled  states  themselves,  and  this  entanglement  has 
interesting  consequences  if  the  cavities  are  fairly  highly 
excited,  as  Gea-Banacloche  [10]  and  Phoenix  and  Knight 
[11]  originally  pointed  out.  Here  we  are  interested  in  non¬ 
local  atom-atom  entanglement  between  the  two  spatially 
separated  sites  on  a  qubit  lattice,  as  suggested  in  Fig.  2. 

In  the  general  context  of  entanglement  we  note  that 
there  is  no  accepted  and  practically  workable  criterion  for 
determining  separability  of  arbitrary  four-particle  states. 
Our  purposes  will  be  satisfied  by  working  with  the  two- 
qubit  atomic  states  obtained  from  the  time-evolving  four- 
qubit  state.  All  familiar  measures  agree  about  separabil¬ 
ity  in  the  two-qubit  domain  of  entanglement.  That  is, 
entropy  of  formation,  Schmidt  number,  tangle,  negativ¬ 
ity,  and  concurrence  are  numerically  somewhat  different, 
but  in  the  two-qubit  domains  of  their  applicability  they 
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are  in  full  agreement  when  they  signal  complete  separa¬ 
bility,  i.e. ,  the  lack  of  entanglement. 

We  adopt  Wootters’  concurrence  [12]  as  our  measure 
in  this  discussion,  mainly  for  its  convenient  normaliza¬ 
tion:  1  >  C  >  0,  where  C  =  0  indicates  separability 
(zero  entanglement)  and  C  =  1  means  maximal  pure 
state  entanglement,  as  in  a  Bell  state;  and  simplicity  of 
calculation: 

=  (io) 

where  the  quantities  A 8-  are  the  eigenvalues  in  decreasing 
order  of  the  auxiliary  matrix 


It  is  easy  to  see  that  a  =  ±7r/4  reproduces  the  two  pure 
Bell  states. 

The  initial  state  for  the  atoms  plus  cavities  is  therefore: 

|$(0)>  =  |0a,06>  (15) 

=  {cosa\eA,eB)  +  sina| gA,gB))  (g>  |0a,  0&). 

In  terms  of  the  dressed  eigenstates  given  above  (5),  we 
can  rewrite: 

|ezi,0a)  =  c\ipf)  -  s|Vq) 

\9aAo)  =  s\^+)  +  c\^)  and 

\9A,$a)  =  \ipo)-  (16) 


C  =  P{<Ty  ®  <Ty)p*{<Ty  (g>  <Ty),  (H) 

where  p*  denotes  the  complex  conjugation  of  p  in  the 
standard  basis  and  ay  is  the  Pauli  matrix  expressed  in 
the  standard  basis. 

One  finds  that  in  reduction  to  two-qubit  form,  by  trac¬ 
ing  over  the  two  field-mode  qubits,  the  resulting  two- 
qubit  mixed  state  always  has  the  X  form  [13],  where  only 
the  diagonal  and  anti-diagonal  elements  are  not  zero: 

'a  0  0  w 

0  b  z  0 

P~  0  z*  c  0 

_w*  0  0  d 

where  a  +  b+  c  +  d=  1.  The  concurrence  of  this  mixed 
state  is  easily  found  to  be 

0  =  2  max{0,  \z\  —  y/ad,  |w|  —  y/bc}  =  2  max{0,  Q},  (13) 

so  it  is  clear  that  Q,  which  is  the  larger  of  \z\  —  y/ad 
and  | in  |  —  y/bc,  will  be  an  important  quantity.  For  exam¬ 
ple,  Q(t)  obeys  certain  conservation  relations  (whereas 
O  does  not)  in  some  special  cases  because  it  can  be  neg¬ 
ative.  Furthermore,  one  can  see  [14]  that  while  Q  <  0 
implies  a  separable  state  (O  =  0),  the  slightly  stronger 
condition  Q  <  0  implies  both  0  =  0  and  also  that  the 
state  is  mixed  rather  than  pure,  information  not  available 
from  O. 

In  the  two-site  situation  under  consideration  there  are, 
in  principle,  six  different  concurrences  that  can  provide 
information  about  the  bipartite  entanglements  that  may 
arise.  With  an  obvious  notation  we  can  denote  these  as 
CAB,  Cab,  CAa,  CBb,  CAb,  CBa.  Except  for  CAa  and 
CBb ,  they  measure  remote  (non-local)  entanglements. 

IV.  PARTIALLY  ENTANGLED  BELL  STATES 

Aspects  of  the  entanglement  dynamics  of  qubits  in 
the  vacuum  two-cavity  case  have  been  reported  already 
[15-17].  For  background  we  summarize  results  when 
the  initial  states  are  superpositions  of  the  Bell  states: 
|$±))~  |  eA,eB)  ±\9a,9b),  which  we  write: 


(12) 


Thus  the  initial  atom-atom  entangled  state  has  the  form 

cos  a\eA,  0a)  (g>  \eB,tib)  +  sin  a\gA,  0a)  (g>  \gB,0b) 

=  cosa(c\ip^)A  -  s \4>i)A)  <X>  (c\4>i)b  ~  s\'4’i)b) 

+  sinalV’0)/!  ®  |V’°)b-  (17) 

Evolution  in  time  is  easily  followed  through  the  evolu¬ 
tion  of  the  dressed  states: 

W±(t))  =  e~iX±t  1^(0)).  (18) 

Note  that  since  the  combination  of  coefficients  in  |$(0)) 
uniquely  associates  c  with  \i>+),  and  s  with  the  time 

evolution  can  be  transferred  to  the  c  and  s  symbols.  We 
will  henceforth  consider  them  carrying  the  time-evolution 
exponents.  We  will  use  the  notation  c o  =  c(0)  and  s o  = 
s(0)  to  refer  to  their  values  at  t  =  0  (no  relation  to  the 
n  =  0  subscripts  in  Eq.  (8)),  so  that, 

c  =  c(t)  =  coe-;iA+t, 

s  =  s(t)  =  s0e“iA-t,  (19) 

where  A+  and  A-  are  obtained  by  inserting  n  =  1  into 
Eq.  (3).  Then  we  can  write  (temporarily  again  indicating 
explicit  time  dependences  for  the  c’s  and  s’s): 

|$M)  =  COS  a  (c(t)\ipA) A  -  s(t)IV’r)zt) 

®  -  s(t)iV’r)s) 

+  sin  a\ip°) A  ®\ip°)B,  (20) 

where  \^)  will  continue  to  refer  to  the  JC  states  at  t  =  0. 

Now  we  revert  to  the  bare  basis  states  in  preparation 
for  the  tracing  needed  to  calculate  QAB  and  CAB: 

|$(t))  =  cosa^c(t)(c0|ey4, 0a)  +  s0\gA,  la)) 

-  s(t)(-s0\eA,  0a  +  c0\gA,  la)))  <g>  (c(f)(c0|eB,  06) 

+  so\gB,  U))  -  s(I)(-so |e_e,  06  +  c0\gB,  U))) 

+  sin  a\gA,  0a)  (g)  \gs,0b) 

=  cos  a  ^(cc0  +  sso)|eA,  0a)  +  (cs0  -  sc0)|^,  la)^ 

(g>  ((cc0  +  sso)|e_B,  Oft)  +  (cs0  -  sc0)\gB,  U)) 

+  sin  a\gA,  0a)  (g)  \gB,0b)- 


\®ab)  =  cos  a\eA,eB)  +  sina\gA,gB).  (14) 


(21) 


4 


To  get  the  two-qubit  mixed  state  needed  for  calculation 
of  QAB  the  projections  that  are  needed  are: 

(0a,  Of, !$(£))  =  cqsq-(cc0  +  sso)2|e0,  eB)  T  sma\gA,gB) 
(la,  Ofc|$(<))  =  cos  a(cs0  -  sc0)(cc0  +  ss0)\g  a,  eB) 

(0a,  l6|<h(f))  =  cos  a(cc0  +  ss0){cs0  -  sc0)\e  a,  gs) 

<l«,l&l*(t)>  =  0.  (22) 

These  show  that  the  AB  mixed  state  has  the  X  form 
mentioned  in  the  previous  Section: 


'a  0  0  ic' 

0  6  0  0 
0  0  e  0  ’ 

.  w*  0  0  d  . 


(23) 


for  which  the  concurrence  has  the  stated  form 

CAB  =  2  ?»aa:{Q,  |«i|  —  \/bc}.  (24) 

and  we  easily  find  the  following 

| w |  =  |  sin  q-cqs  o-|(Cq  +  Sq  +  2cqSq  cosJt), 

6  =  c  =  cos2  a|cco  +  sso  |2  |cso  —  sco|2 
=  COS2  «(Cq  +  Sq  +  2CqSq  cos  St) 
x  CqSo(2  —  2cos St ),  (25) 


where  S  =  \J A2  +  4 g"2 . 

For  simplicity  we  will  evaluate  this  in  the  resonance 
case,  6n  =  7t/2,  where  co  =  so  =  l/y2  and  A  =  0.  Then 
we  find 


|«)|  —  Vbc  =  —  cos2  cv[2  +  2  cos(gt )][|  tana]  —  sin2(<?t)], 

'  (26) 
from  which  the  expression  for  concurrence  is  found  to  be: 


CAB  '1  rtiu.x{n.  (J  W!|.  (27) 

where  QAB  =  cos2  a  cos2 (gt.  )[|  tana]  —  sin2  (</£)]. 

Fig.  3  shows  that  for  a  yf  7t/4  the  CAB  curves  have  the 
“sudden  death”  feature  [18].  That  is,  the  entanglement 
non-smoothly  becomes  zero  and  stays  zero  for  a  finite 
interval  of  time.  It  is  also  clear  that  the  AB  and  ab  en¬ 
tanglements  are  opposites  in  the  sense  that  one  increases 
as  the  other  decreases  (except  in  the  regions  where  both 
are  zero).  This  has  been  analyzed  in  terms  of  information 
transport,  and  time-invariant  relations  have  been  identi¬ 
fied  [15]  that  connect  all  six  Q(t)  values.  Because  the 
evolution  in  the  ideal  situation  treated  here  is  lossless, 
and  only  finitely  many  states  are  involved,  the  periodic 
nature  of  the  curves  is  not  surprising.  However,  the  sit¬ 
uation  is  more  realistic  when  infinitely  many  final  states 
are  involved  (quantum  open  system  case).  Results  are 
available  in  cases  involving  simple  harmonic  oscillators 
that  interact  with  reservoirs  [19,  20].  In  the  remaining 
sections  we  will  employ  the  techniques  presented  here  to 
extend  the  vacuum  results  into  the  open-system  domain 
by  allowing  the  atoms  to  interact  with  (be  controlled  by) 
more  realistic  fields,  namely  fields  described  by  quantum 
coherent  states. 


FIG.  3:  Time  dependences  of  AB  and  ab  concurrences  for 
three  values  of  the  superposition  angle  a.  Note  that  almost 
all  curves  have  some  time  interval  over  which  C=0,  i.e.,  during 
which  the  underlying  state  must  be  separable.  From  ref.  [15]. 


V.  OPEN-SYSTEM  TWO-QUBIT  THEORY 


Realistic  quantum  control  almost  necessarily  implies 
engagement  of  continuous  variables  and  the  interaction 
of  qubits  or  other  systems  having  a  finite  number  of  states 
with  one  or  more  “large”  systems  with  continuous  quan¬ 
tum  states.  Novel  bipartite  aspects  of  open-system  joint 
qubit  states  have  been  noted  recently  [21].  We  can  use 
the  JC  formalism  described  above  to  enter  this  domain 
by  introducing  coherent  state  fields  at  the  AB  network 
sites.  The  Jaynes- Cummings  interaction  remains  rele¬ 
vant,  and  we  again  focus  on  pure  storage  memory  nets, 
with  sites  entangled  but  not  interacting,  as  in  Fig.  2. 

We  retain  almost  all  of  the  simplifying  approximations 
made  earlier,  and  add  one  more  by  taking  the  two  fields, 
now  modeled  as  coherent  state  fields,  to  have  the  same 
average  photon  number  n.  However,  one  important  sim¬ 
plification  that  we  relied  on  previously  must  now  be  dis¬ 
carded.  The  coherent-state  fields  have  many  occupied 
photon-number  states,  so  the  cavities  will  not  be  qubits. 
We  assume  initial  atom  entanglement  in  the  form  of  a 
pure  Bell  State: 

I*ab(0)>  =  (\eg)+\gt))/V2,  (28) 

and  write  the  field  state  as  the  coherent  state  product 
|tta6(0))  =  |a>O|a).  (29) 

As  a  result  our  initial  state  for  the  whole  system  is, 


|tfta*(Q)>=  |*AB(0))O|a)o|a>.  (30) 


The  coherent  states  are  given  by  the  familiar  Fock  state 
expansion 


]®>  =  H 

n= 0  n=0 


\n). 


(31) 


Because  a  coherent  state  provides  an  infinite  num¬ 
ber  of  Fock  states,  the  photonic  density  matrix  be¬ 
comes  infinite-dimensional,  and  the  joint  AB  dynamics 
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extremely  complicated.  Fully  numerical  analysis  is  pos¬ 
sible,  but  the  insights  from  analytic  results  are  highly 
desirable.  We  have  found  a  key  step  permitting  ana¬ 
lytic  progress.  This  is  an  apparently  drastic  simplifica¬ 
tion  of  the  continuous  state  spaces  of  the  two  held  modes. 
We  assume  that  it  is  satisfactory  to  replace  |a)  by  |h). 
This  Ansatz  is  at  least  weakly  supported  by  the  knowl¬ 
edge  that  photon  number  in  a  coherent  state  is  Poisson- 
distributed  and  relatively  tightly  centered  around  n  when 
n  1.  Thus  we  represent  the  initial  held  state  by  the 
single  product  Fock  state  |h)  (g)  |h).  Note  that  during  the 
JC  interaction  the  photon  number  in  either  held  mode  a 
or  b  can  then  be  n  or  n  +  1  or  n  —  1.  Thus  the  held  modes 
are  no  longer  convenient  qubits  as  they  were  in  Sec.  IV. 

Under  the  Ansatz  mentioned,  by  tracing  the  held  mode 
states,  we  hnd  that  the  reduced  density  matrix  for  the 
qubits  becomes: 


where 

Kmn  —  iCn  + 1  \e ,  6 ,  71 ,  TTl  1)V  Cn  + 1  1 6 ,  g ,  71 ,  771) 

-  Nn+iNnl#,  e;  n  +  1,  m)  -  iSn+iC„ \g,  g;  n+l,m) 

-  iSnCm+i  |e,e;n  —  l,m+  1) 

-  S„Sm+i\e,g;n  -  l,m  +  1) 

+  C„Cm+i\g,e;n,m+  1) 

-  iC„Sm+i\g,g;n,m+l),  (38) 

where  Cn  =  cos(gt^/n)  and  Sn  =  sin(gt^/n).  The  density 
matrix  for  the  system  is  then  given  by 

Ptot  =  |tftot(*))(|tftot(*)|.  (39) 

and  the  reduced  density  matrix  for  the  atoms,  Pab,  is 
given  by, 

PAB  —  Tr)n!m)  Ptot  •  (40) 

Having  used  the  Fock  state  shortcut  to  obtain  (32), 
we  avoid  using  it  further  now  and  calculate  the  elements 
z,  a,  d  of  pab  for  the  coherent  state.  The  z  term  is  given 
by  the  doubly  inhnite  summation, 


where  we  have  used  the  standard  two-qubit  basis 
[ee,  eg,  ge,  gg\.  The  elements  indicated  by  x  are  zero 
because  of  the  equal-h  simplihcation.  Thus,  under  the 
assumptions  mentioned,  p  is  again  of  A-type  [13]). 

For  this  A-type  p,  concurrence  turns  into: 

C(p)  =2  max[  0,  \z\  —  \fad  ].  (33) 


A^A^nCnCn  +  iCrnCrn  +  i 


AnAn  —  i  Am  Am_j_i  Sn  CfiAiCmSmAi 
T  An  An_2  Am  Am_|_2  Sn  Sn  —  iNm-j-i  SmA2 
AnAn  —  iAmAm_|_i  Sn  Cn  —  iNm-j-i  Cm-\-2 


}.  (41) 


Similarly  the  series  summations  for  a  and  d  are: 


The  coherent  fields  induce  time-dependent  change  of  the 
elements  a  and  d,  and  their  growth  and  any  decline  of  z 
will  cause  entanglement  to  decrease. 

In  order  to  calculate  the  time  evolution  of  this  state 
under  the  JC  Hamiltonian  we  need  to  calculate  the  time 
evolution  of  the  the  states  of  the  individual  sites,  i.e. , 
|e)  (g>  |n)  and  |g)  (g)  |n)  for  all  n.  The  time  evolution  of 
these  states  for  site  A  (and  similarly  for  site  B)  is  given 

by, 

elHlt\ e;n)  =  cos(gt\/ n  +  l)|e;  n) 

—  ism(gt\/n  +  l)|g;  n  +  1)  (34) 

=  cos{gt^)\g-,n) 

—  i  sm(gt^/n)  |e;  n  —  1).  (35) 


The  inhnite  extent  of  these  summations  of  course  reflects 
the  fact  that  we  are  dealing  with  a  quantum  open  system, 
by  having  coupled  the  qubits  to  an  inhnite  state  space. 

The  sums  above  cannot  be  completed,  but  excellent 
analytic  approximations  can  be  found  for  coherent  states 
that  are  even  only  moderately  strong,  i.e.,  a  >  10.  We 
will  use  the  familiar  Stirling  formula  for  n\, 


n!  =  \/2tt nnne  n , 


(44) 
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and  Euler’s  formula  to  approximate  the  terms  in  the  sum¬ 
mations  above  by  integrals.  We  begin  by  approximating 
the  terms  like  A„An+iAmAm-i  with  A^A^,  which  in¬ 
troduces  an  error  of  order  1  jn  near  the  Poisson  peaks 
n  «  m  «  n.  Then  we  get  for  z 

Z  ~  \{{Y.XnCnCn  +  l)2 

n 

-  2{j2A2nA2mSnCn+1CmSm+1 ) 

n,m 

+  (^^5n5n+1)2}.  (45) 

n 

In  the  same  way,  Eq.(42)  and  Eq.(43)  become, 
a  -  (C£AnC^  (j2An'Sn)  +  (E AlSnCny,  (46) 

n  n  n 

d=  (j2A^Cn){j2AnSn)  +  '  (47) 

n  n  n 

Note  that  Eq.  (46)  and  Eq.  (47)  imply  that  a  =  d 
whenever  our  large  n  approximation  is  valid.  Now  we 
rewrite  C'nCn+i  as, 

CnCn+ 1  =  ^  jcos  [gt(y/n  +  \/n  +  1)] 

+  COS  \gt(^/n  +  1  -  y/n)\  },  (48) 

and  use  the  peaked  nature  of  An  to  focus  on  those  terms 
near  to  n  to  introduce  the  approximation 


Now,  using  the  identities 

Cl  = 


2  1  +  cos(2gtyd) 


_  1  COS (2 gty/n) 


SI  = 


we  can  rewrite  a  and  d  as 
1 
4 


=  d-\\l~{^Al  cos(2gty/n)J 

n 

+  (J2A™  sin(2 9tVn)] 


Then  Eqs.  (52)  and  (54)  lead  to: 


+2(E4;  cos (2 gty/n)^ 

n 

“KXE"  siii(2^-Cn))  -1 

n 


(53) 


(54) 


(55) 


We  can  calculate  the  sums  involved  here  by  rewriting 
them  as  integrals,  treating  the  integer  n  as  continuous, 
again  relying  on  the  large  value  of  n.  The  first  two  inte¬ 
grals  we  need  to  calculate  are, 


h  = 


/  A2  cos(  d  )dn,  and 

Jo 


(56) 


Vn  1  =  \fn  + 


(49) 


h 


o 


A2n  sin( 


(57) 


which  compresses  CnCn+i  to 
1 


CnCn-\. i  —  ^  |  cos 


(2 gty/n)  +  cos  (^=)  .  (50) 


Similarly, 


5n5n+i  =  i[cos  -  cos{2gt^/n) 

l  L  \2^/n/ 


S„cn+ 1  =  i  [sin(2^Cn)  -  sin  (^=^) 

S'n+iCn  =  i  [sin (2gt y/n)  +  sin  {^=) 


With  these  results  we  can  simplify  z  further: 


(51) 


n  v 

+  (^2Al  cos(2 gty/n) 

n 

-{J2A™  sin(2 9tVn) 


^2cos(  gt 


E  A 2  sii 


,_gi_. 

'2  y/n' 


(52) 


We  will  combine  these  integrals,  Ii  +  il2  =  1 12  in  order 
to  work  with  the  exponential  of  igt/2*Jn.  This,  together 
with  Stirling’s  approximation  on  A2,  and  the  abbrevia¬ 
tion  t  =  gt,  leads  to 

r°°  rv^n  pn 

I12  =  /  e~a  ^=—eiTC^dn.  (58) 

Jo  y/2^nn  V  7 

The  saddle  point  method  is  appropriate  to  calculate 
this  integral,  and  some  details  are  reserved  for  the  Ap¬ 
pendix.  The  expression  for  Ii  +  il2  is  found  to  be, 


h2  =  exp 


32a4 


e ir/2a 


(59) 


Then  helpful  cancellations  can  be  identified,  and  we 
obtain  an  approximate  expression  for  \z\  —  y/ad: 


\z\-Vlld  =  I[e-C‘2/16n2  _ 

+  ^[E^2  cos(2^c«)] 

n 

~  7}[J2A™  sin(2^C«)]  •  (60) 
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The  summations  in  (60)  involving  cos(2 gt^fn)  and 
sin  (2 gt\/n)  can  be  combined  into  a  single  exponent  con¬ 
taining  the  argument  2 igty/n,  which  is  similar  to  that  in 
(58),  but  the  resulting  saddle  point  analysis  is  more  com¬ 
plicated  because  now  v/ n  is  in  the  numerator.  Details  are 
relegated  to  the  Appendix,  with  the  working  result: 


z  — 


-1  +  e 


-t2/ 2 


y  JL 

2irk 


exp 


2(r  —  27tA>Q')2  \ 
1+7 r2&2  / 


x  cos[4o-(r  —  'Ijrka)] 


(61) 


In  writing  this  last  equation  we  have  used  the  fact  that 
around  r  =  2i rka  only  the  term  with  the  corresponding 
k  gives  a  significant  contribution  to  the  sums.  The  con¬ 
tribution  to  t  =  2irka  from  any  other  k'  is  proportional 
to  exp  ^  —  47r2Q'2(k  —  fc')2/[l  +  7r2(k')2]J,  so  it  decays 
exponentially  with  the  distance  from  k. 


VI.  OVERVIEW  AND  IMPLICATIONS 

We  have  demonstrated  the  utility  of  a  novel  short-cut 
approach  to  determine  the  open-system  evolution  of  en¬ 
tanglement  under  coherent-state  interaction  with  two  re¬ 
mote  qubits.  The  A’-state  simplification  and  Fock-state 
Ansatz  we  introduced  work  together  in  a  surprisingly  ac¬ 
curate  way.  Some  details  of  the  calculations  themselves 
have  been  included  for  inspection,  and  we  can  summa¬ 
rize  them  by  greatly  over-simplifying  the  multiple  stages 
involved  as  follows.  The  evaluation  of  concurrence  for 
n  1  here  is  generically  the  same  as  that  presented 
for  zero  detuning  qubit  inversion  in  the  original  discus¬ 
sion  of  quantum  revivals  [23] ,  so  one  should  be  prepared 
to  find  revival  behavior  here,  and  strong  revivals  of  en¬ 
tanglement  are  apparent  in  Fig.  4  below.  They  show 
both  the  predictions  of  the  analytic  expressions  given 
above,  and  also  the  results  of  a  numeric  check  on  the 
Fock  state  Ansatz  adopted  to  initiate  the  analytic  calcu¬ 
lations.  That  is,  Fig.  4  shows  our  analytical  results  for 
the  function  2(  \z \  —  \fad)  in  comparison  with  the  complete 
numerical  evaluation  of  concurrence  without  making  the 
A-type  approximation.  The  agreement  is  not  perfect, 
but  the  discrepancies  are  illuminating. 

In  particular,  the  results  for  entanglement  contain  de¬ 
tails  not  present  in  well-known  results  for  inversion,  and 
these  details  appear  to  work  favorably  for  application  to 
quantum  memory  management,  as  to  both  storage  and 
recovery.  This  is  not  difficult  to  see,  as  we  now  indicate. 

There  is  an  important  scale  issue  to  note  in  these  fig¬ 
ures.  Both  curves  in  Fig.  4  show  long  periods  of  ESD 
interrupted  by  intervals  of  coherent  response,  i.e. ,  recov¬ 
ery  of  entanglement.  While  it  is  not  perfect  recovery  it  is 
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FIG.  4:  The  analytical  and  numerical  results  for  concurrence 
of  two  qubits  exposed  to  two  quantum-coherent  driving  fields 
with  a  =  10.  The  evident  revivals  are  predicted  reasonably 
well  by  the  approximate  formula  (61).  Analytical  results  are 
for  the  A-type  p  while  the  numerical  ones  are  for  the  original 
P- 


surprisingly  robust.  The  envelope  of  the  revivals  decays 
relatively  slowly,  as  given  by  the  formula 

J_  _  1  -exp(-r2/16a4) 

7 rk  2  1  ] 

where  k  is  the  revival  number. 

Returns  of  concurrence  to  zero  (ESD  events)  within  the 
analytic  revivals  are  not  present  in  the  numeric  revivals. 
The  extent  of  this  is  seen  by  inspecting  a.  revival  envelope 
in  detail,  as  is  done  in  Fig.  5.  Close  to  the  revival  regions 
the  analytical  results  show  rapid  oscillations  with  period 
t  =  ir/(2g\/fi)  which  is  half  the  corresponding  period 
in  the  inversion  case.  The  analytic  formula  retains  the 
entanglement  death  and  rebirth  episodes  that  occur  on 
the  rapid  Rabi-period  scale,  as  were  shown  in  Fig.  3 
and  have  been  discussed  in  the  literature  repeatedly  for 
few-photon  excitation.  By  contrast,  the  numeric  results 
show  a  smoothed  version  without  rapidly  recurring  ESD 
events.  Even  the  zero  revival,  i.e.,  the  period  that  is 
referred  to  as  the  Cummings  collapse  in  the  inversion 
literature,  shows  no  ESD  events  within  it. 

One  can  conclude  that,  in  principle  a.t  least,  several 
aspects  of  quantum  memory  management  can  be  fur¬ 
ther  developed  for  remote  qubit  entanglement  using  the 
interactions  employed  here.  For  example,  we  see  that 
quantum  joint  memory  can  be  “hidden”  without  easy  de¬ 
tectability.  In  the  ESD  intervals  the  stored  entanglement 
is  not  evident  either  locally  or  globally.  Moreover,  these 
long  intervals  have  two  important  management  charac¬ 
teristics  -  they  have  closely  predictable  durations,  and 
they  don’t  substantially  deplete  the  degree  of  the  stored 
entanglement. 

Finally,  we  can  comment  on  the  quasi-periodic  mod¬ 
ulations  evident  in  the  numeric  details  in  Fig.  5.  For 
contrast,  we  show  in  Fig.  6  below  the  numeric  revival 
det  ails  for  two  other  values  of  coherent  state  photon  num- 
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FIG.  5:  Details  of  the  first  revival  shown  in  the  preceding 
figure.  The  deep  modulations  in  the  analytic  envelope  are  not 
present  in  the  numerical  envelope,  which  retains  entanglement 
robustly  through  the  revival  event. 
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FIG.  6:  Similar  to  Fig.  5  above,  except  that  here  a  =  5  and 
6,  rather  than  10. 


ber:  n  =  25  and  36.  In  those  graphs  modulations  also 
appear,  but  with  different  main  periods,  their  frequency 
increasing  linearly  as  a  increases,  viz.,  n  main  modula¬ 
tion  periods  for  a  =  n  over  a  unit  interval  in  gt/n.  These 
modulations  would  be  further  smoothed  in  realistic  ap¬ 
plications  where  differences  would  have  to  be  expected 
between  the  n  values  at  the  two  sites  being  managed, 
thus  promoting  rather  than  degrading  the  desirable  fea¬ 
tures  mentioned. 
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VIII.  APPENDIX 


Here,  we  sketch  the  saddle  point  analysis  of  the  re¬ 
maining  sums  in  Eqn.  (60).  They  can  be  combined,  and 
converted  to  an  integral: 

pOO 

h 4=  /  exp  [a2f(n)\d.n,  (63) 

Jo 

where 

f(n)  =  —  (2n .  lncv  —  —  ln(27rn) 
mV  2 

—  11  In!  n  i  +  11  +  2  iTy/n)  —  1.  (64) 

In  the  saddle  point  method,  an  integral  of  type 
fa  eM^^dx,  where  M  is  a  large  number  and  f(x)  is  a 
twice-differentiable  function,  is  approximated  by, 


eMf{x)dx 


2? r 


M\f"(x0)\ 


pMf  (ro 


(65) 


where  xo  is  the  global  maximum  of  f(x).  In  our  appli¬ 
cation  Eq.  (63)  can  be  written  in  this  form  if  we  choose 
M  =  cr2  and  f(n)  as  given  above.  For  the  saddle  point 
maximum  of  f(n)  we  need  to  find  the  point  no  where 

/'(no)  =  0: 


/'(no)  =  —(21110--- - ln(fi0)  +irn0  =0.  (66) 

Assuming  that  |tto|  is  large,  we  can  replace  ln(no)  — 
ln(o2)  +  irn0  1/'“,  and  then  letting  no  =  pet6  gives 

In  p  +  id  =  ln(o2)  +  iTp~1/2e~te/2,  (67) 

where  equating  the  real  and  imaginary  parts  of  both  sides 
gives, 

ln(p)  =  ln(o2)  +  Tp-1!2  sin(0/2) 

6  =  rp~1/2  cos(6>/2).  (68) 

For  r  =  27t fco,  where  k  is  an  integer,  the  equations 
above  become, 


p  =  a2  and  6  =  (  —  l)k2jTk.  (69) 


Now,  let  t  =  To  +  Ar  and  0  =  do  +  AH,  where  To  =  2irka 
and  0 0  =  2irk.  Assuming  that  both  Ar  and  AH  are  small, 
we  can  write 


sin(H/2)  =  (  —  l)fc-^-  and  cos(H/2)  =  (  —  l)fc.  (70) 


Then  Eq.  (68)  turns  into, 


P 

0 


o2[l  +  (  — l)fc7rfcAH] 

wkAd  h  At 

£0-  — +  M)  —  • 


(71) 


'o  - 


2 


9 


In  order  to  arrive  at  the  second  line  we  have  used, 


!-(-!)* 


irkA9y 

2  )’ 


(72) 


Again,  using  9  =  90  +  A 9  and  r  =  To  +  At,  writing  A 9  in 
terms  of  At  and  ignoring  the  terms  after  second  order, 

/m{a2/(rio)}  =  (—l)k  Alirka2  +  Via 


in  the  second  line  of  Eq.  (68)  and  we  have  ignored  the 
terms  with  AtA9. 

Now,  the  last  two  terms  in  the  second  line  of  Eq.  (71) 
are  just  A 9,  so, 


2  a 


Thus,  we  have, 


A  9  = 


(-1)* 


a(l  +  ir2k2) 
Using  this  in  Eq.  (71),  we  obtain 

1 


At, 


p  =  cr[l  +  7T  k 
9  =  90  +  (-!)' 


a(l  +  ir2k2) 

,  1 


Ar] 

At. 


a(l  +  7 x2k2) 

Now,  by  inserting  no  =  pe10  in  Eq.  (89),  we  find 


2f(n o)  =  -peie  ln(p/a2) -ln(p1/2) - 

pei0(  1  -  i9)  +  liTP1/2eie/2 
—  a2  —  ln(V27r). 


i9 


(73) 


(74) 


(75) 


(76) 


The  real  part  of  this  equation  is  given  by, 

Re{a2f(no)}  =  —  p  cos  9  ln(p/a2)  —  ln^1/2)  +  p9  sin  9 
+p  cos  9  —  2 rp1!2  sin((?/2) 

-a2 -ln(V2^).  (77) 


Using  9  =  9 o  +  A 9  and  retaining  only  the  terms  upto  the 
second  order  in  A  9, 


Re{a2f(n0)}  =  -p[  1  -  In  (y) 

-  In (p1/2)  +  pA9(90  +  A 9)  +  p  [l  -  ] 

—  (— l)kTp1!2A9  —  a2  —  ln(V27r).  (78) 

Writing  A 9  in  terms  of  Ar  and  ignoring  the  terms  after 
second  order, 


fle{“3/(""))  - -irW(lrfi)(iT)2 

—  In)/?1/2)  —  In)/!?1/2)  —  ln^1/2)  —  ln(V27r).  (79) 

The  imaginary  part  of  a2 /(no)  is, 

Im{a2 /(no)}  =  —  sin  9p\n(p/a2)  —  9/2 
—p9  cos  9  +  psin  9  +  2 rp1/2  cos(0/2).  (80) 


1 


Ar 


2a(l  +  ir2k2) 
irk  f  3 


1  +  7 x2k2  L 2 ( 1  +  ir2k2) 


( A7 


(81) 


For  k  =  0  the  (Ar)2  part  of  this  equation  vanishes.  For 
k  =  1,2,...  this  part  can  be  ignored  as  well.  Thus  we  are 
left  with, 

Im{a2  f(no)}  =  (— l)k(lirka2  +  2aAr).  (82) 

We  also  need  to  calculate  \/2ir  /  a2\f"  (no)\.  By  using 
Eq.  (66),  we  find 


,,,  ,  1/1  1  irn0 


-3/2 


a2  \2 n2  no 
which  has  two  forms:  for  r  =  0, 


2t r 


Y  «2|/"(«o)| 
and  for  r  =  27T ka  (k=  1,2,...): 


=  V  27TO, 


2? r 


27 ra 
irk 


V  «2l/"<>o)|  ’ 

Then  the  integral  in  Eq.  (63)  is  given  by, 


Um  — 


2tt 


a2|/"(«o)| 


oa2/(n0 


(83) 


(84) 


(85) 


(86) 


In  order  to  find  the  value  of  this  integral  we  should  add 
the  contributions  from  all  k  =  0,  1,  2,  ....  As  a  result,  the 
integral  is, 


U4  =  e~T2/2e2io 


I  1  r  (r  —  2irka)2 

£  V^expr 

fc  =  l,2,... 

cos[2a(r  —  2irka)]. 


1  +  ir2k2 


(87) 


Now,  inserting  the  values  of  all  four  integrals  into  Eq. 
(55), 

z  —  \fad  =  —  |  exp  ^  —  yy — —  1  +  e~T  cos(4ar) 

2  1  r  /  2(r  —  2irka)2 

Z  I  exp 


7T  A — '  k 

k  =  1,2 


1  +  ir2k2 

[4a(r  —  27rfca)]J  j-.  (88) 


Writing  this  last  equation  we  have  used  the  fact  that 
around  r  =  2irka  only  the  term  with  the  corresponding 
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k  gives  a  significant  contribution  to  the  squares  of  the 
sums  in  Eq.  (55).  The  contribution  to  r  =  2n ka  from 
any  other  k'  is  proportional  to  exp(— 4ir2a2(k  —  k')2 /[  1  + 
7r2(fc')2]),  so  it  decays  exponentially  with  the  distance 
from  the  peaks  identified  with  integer  k.  Thus,  while 
taking  the  squares  in  Eq.  (55)  we  can  ignore  the  cross¬ 
terms. 

Now,  we  turn  to  the  integral  given  in  Eq.  (58).  We 
will  continue  to  use  the  saddle  point  method.  This  time 
we  need  to  find  the  maximum  of  the  function: 

f(n)  =  —  (2nln(a) - ln(27rn) 

az  2 

IT 

-nln(n)  +  n  +  — -=)  -  1,  (89) 

2  V« 

For  this,  we  need  to  find  the  point  no  where  f  (no)  =  0: 
/>o)  =  2^_ln(n°)_Tn“3/2')  =  0> 


Writing  p  in  terms  of  9  and  retaining  only  the  terms  up 
to  82 , 


a2  f(n  o)  =  —ln(a2)  —  ln(a)  —  ln(V27r) 
36»2  a292  t0 


4  2  4a 

(  T  6  7 T02 


so  that  by  inserting 


2  16a/’. 

for  6,  we  obtain 


(97) 


a2f(n0) 


ln(a2)  —  ln(a)  —  ln(V27r) 
r2  3r2 
32a4  16a6 

8^  ~  256^)’  (98) 


Again  assuming  that  |rio|  >>  1, 

In (n0)  =  ln(a2)  -  yrig  ' 3/2 


(91) 


In  order  to  find  the  final  form  of  Eq.  (65)  we  need  to 
calculate  |/"(rio)|  as  well. 


Letting 


i» 

n0  =  pe  , 


then, 


(92) 


/>o) 


1/1  1  3ir  —5/2 

^\2ni  ~~0+~n° 


(99) 


it  S/9  S/9 

In (p)+id  =  ln(a2)  -  ^-p“3/2[cos(y)  -  isin(y)],  (93) 

and  by  matching  the  real  and  imaginary  parts  of  both 
sides,  we  find  two  coupled  transcendental  equations: 


Retaining  the  terms  up  to  a  2  in  the  parentheses,  we 
have  |/" (no) |  —  /pr-  Then  we  use  Eq.  (98)  and  Eq.  (65) 
and  insert  a2  for  M,  to  get 


ln(p)  =  ln(a2)-y  3/2sin(y), 

9  =  ~^P~3/2  cos(y).  (94) 

We  are  going  to  retain  only  the  terms  up  to  second  order 
in  rp-3/2, 


2tt 


a2|/"(n0)| 


ocx2f(n0 


exp 


3r2  n 
32a4  16a6 / 


ir 

8a3 


7  ir3 
256a 


For  r  ~  a2  we  can  retain  the  first  terms  in  the  parenthe¬ 
ses  and  ignore  the  rest  since  a2  >>  1.  Thus,  we  finally 
find 

With  this  restriction,  and  after  inserting  n  =  pe10  into 
Eq.  (89),  we  find 


T 

4a3 


S/92 

and  p  =  a2(l  4 - ).  (95) 


a2/(n0) 


—ln(a2)  —  ln(a)  —  \h(\/2tt) 


4 

ir 


i0  t 
-2+?6 


p-1l2e~i0l2  —  a2. 


Ii2  =  exp 


32a4 


e%r  /  2a 


(101) 


(96) 
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I.  INTRODUCTION 

A  quantum  system  of  interest  may  be  large  or  small, 
but  its  background,  called  an  environment,  is  almost 
always  complex,  and  is  often  represented  by  a  bath  of 
bosons  or  fermions,  or  by  classical  random  fields.  In  all 
these  cases,  system  dynamics  is  described  by  a  quantum 
master  equation  that  governs  the  evolution  of  the  reduced 
density  matrix  of  the  system.  In  the  current  decade,  de¬ 
coherence  dynamics  of  entangled  quantum  systems  un¬ 
der  the  influence  of  environmental  noises  has  been  ex¬ 
tensively  discussed  in  different  contexts  involving  atoms, 
ions,  photons,  quantum  dots,  and  Josephson  junctions, 
to  name  several.  This  is  all  related  to  new  regimes  of 
information  processing,  such  as  quantum  cryptography 
and  quantum  computation  [1],  An  important  category  of 
such  research  has  treated  the  fascinating  domain  where 
entanglement  of  qubits  evolves  even  though  the  qubits  do 
not  interact,  even  indirectly.  An  example  is  sketched  in 
Fig.  1  and  we  restrict  our  attention  here  to  this  category. 

In  experimental  environments  an  entangled  system 
may  be  exposed  to  vacuum  noise,  phase  noise,  thermal 
noise,  and  various  classical  noises,  as  well  as  mixed  com¬ 
binations  of  noises.  A  number  of  idealized  models  have 
provided  new  insights  by  allowing  entanglement  evolu¬ 
tion  to  be  followed  by  solving  the  appropriate  quantum 
master  equation  (see  Zyczkowski,  et  al.  [2],  Da.ffer,  et  al. 
[3],  as  well  as  [4-9]).  Most  research  on  entanglement  dy¬ 
namics  has  been  focused  on  ambient  noises  from  environ¬ 
ments  that  obey  the  Markov  (no  memory)  assumption. 
Recently,  there  is  growing  interest  in  the  non-Markovian 
entanglement  dynamics  for  both  discrete  and  continuous 
quantum  systems  (see  [10-15]  and  the  overview  in  [9]). 


In  truth  every  environment  is  non-Markovian.  Non- 
Markovian  noise  was  a  repeated  theme  in  the  research  of 
Krzysztof  Wodkiewicz  [16-27],  and  we  present  our  find¬ 
ings  as  a  contribution  to  his  scientific  memory. 

As  far  as  we  know  there  are  no  fully  systematic  in¬ 
vestigations  of  non-Markovian  noises  or  of  their  effect 
on  the  coherence  dynamics  of  non-interacting  spin  sys¬ 
tems.  In  particular,  a  perturbative  theory  leading  to  the 
Markov  approximation  is  still  lacking.  The  purpose  of 
this  paper  is  to  present  a  study  of  such  problems  in  the 
simplest  form.  We  will  consider  classical  non-Markovian 
noises,  modelling  them  as  so-called  Ornstein-Uhlenbeck 
processes,  and  derive  the  consequences  for  entanglement 
dynamics.  This  can  be  considered  an  extension  of  our 
earlier  note  on  entanglement  sudden  death  (ESD)  under 
classical  Markov  noises  [28]. 


II.  THE  KUBO-ANDERSON  MODEL 
EXTENDED  TO  TWO  QUBITS 

We  consider  an  entangled  pair  of  spins  both  of  which 
are  subject  to  frequency  fluctuations  that  are  random 
[29].  We  adopt  a  model  for  these  fluctuations  that 
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FIG.  1:  Sketch  of  remote  qubits  in  a  quantum  memory  net, 
where  the  dashed  lines  indicate  entanglement,  but  not  inter¬ 
action. 
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treats  them  as  caused  by  noisy  environments  described  by 
Ornstein-Uhlenbeck  processes.  This  well-known  Gaus¬ 
sian  noise  model  is  non-Markovian  in  the  general  case  but 
has  a  well-defined  Markov  limit.  To  focus  exclusively  on 
the  effects  on  the  entanglement  of  the  spins  as  it  arises 
from  the  noise,  we  assume  the  spins  to  be  affected  sepa¬ 
rately  by  separate  environments,  and  not  to  interact  with 
each  other  in  any  way,  especially  not  through  the  noises. 
Thus  the  spins  could  be,  for  example,  remote  compo¬ 
nents  of  a  quantum  memory  net  (as  in  Fig.  1)  under 
steady  attack  by  weak  local  noise.  This  compromises  the 
preservation  of  their  entanglement. 

The  Flamiltonian  of  the  two-spin  system  can  then  be 
formally  written  as  (we  set  h  =  1): 


B{t ) 

- - - <T, 


(1) 


where  £2^(2)  and  £2^(2)  are  the  independent  fluctuations 
of  the  spin  transition  frequencies  (level  spacings).  They 
have  the  mean  value  properties 


M[£2,-(2)] 

M[£2i(2)£2.-(s)] 


0,  (2) 

a(t  —  s) 

i  =  A,B,  (3) 


where  M[-\  stands  for  the  statistical  mean  over  the  noises 
£2^(2)  and  £2^(2).  Note  that  7  is  the  noise  bandwidth, 
and  7-1  =  rc  defines  the  environment’s  finite  correlation 
time  of  the  noise.  For  simplicity,  we  will  take  the  noise 
properties  to  be  the  same  for  A  and  B  (e.g.,  =  T#  = 

T),  although  independent.  In  the  limit  7  — >  00,  Ornstein- 
Uhlenbeck  noise  reduces  to  the  well-known  Markov  case: 


a(t,  s)  =  TS(t  —  s).  (4) 

For  the  total  system  described  by  the  Hamiltonian  (1), 
the  stochastic  Schrodinger  equation  is  given  by 

(5) 

The  explicit  solution  for  the  stochastic  Schrodinger  equa¬ 
tion  can  be  readily  obtained  in  terms  of  a  stochastic  uni¬ 
tary  operator: 


m))  =  u(t,nA,nB)\m),  (6) 

where  the  stochastic  propagator  U  ( 2 ,  £2^4,  £2#)  is  given  by 

lJ(t,ClA,ClB)  =e-^o  (7) 

The  reduced  density  matrix  for  spins  A  and  B  is  then 

obtained  from  the  statistical  mean 

p(t)  =  Mm)mm-  w 


readily  derived  from  the  stochastic  Schrodinger  equation 


[30]: 

dp 

dt 

|'*(2p  ofpof  <rfp<rf), 

(9) 

where 

G(t)  = 

J  a(t  —  s)ds  =  ^(1  —  e_7t), 

(10) 

The  memory  information  of  the  environmental  noises  is 
encoded  in  the  time-dependent  coefficient  G(t)  where 
rc  =  I/7  characterizes  the  environmental  memory  time. 
In  the  Markov  limit  rc  — y  0  (7  — >  00),  when  G(t)  — >■  T /2, 
equation  (9)  reduces  to  the  well-known  Markov  master 
equation  in  the  presence  of  dephasing  noises. 


III.  EXACT  SOLUTIONS  FOR  QUANTUM 
EVOLUTION 

Solutions  of  master  equations  for  the  noisy  evolution  of 
two-spin  density  matrices  in  terms  of  the  Kraus  operator- 
sum-representation  have  been  given  before  (see,  for  ex¬ 
ample,  [28,  31,  32]).  In  many  cases  of  physical  interest, 
the  Kraus  representation  allows  a  transparent  analysis  of 
entanglement  decoherence  without  invoking  the  explicit 
forms  of  the  initial  conditions.  In  what  follows,  we  will 
use  the  fact  that  for  any  two-spin  initial  state  p  (pure  or 
mixed) ,  the  evolution  of  the  reduced  density  matrix  can 
be  written  compactly  as 

P(*)  =  (n) 

where  the  Kraus  operators  satisfy  ^2  =  1  for 

all  2. 

In  order  to  derive  the  desired  Kraus  operators  for  the 
reduced  density  matrix  we  begin  by  noting  that  the  so¬ 
lution  for  just  spin  A  can  be  written: 

\iP(t))  =  u(nA,t)\ip(o))  (12) 

where 

U (£2^,1)  =  exp  [—iF(t)crz]  (13) 

with  the  stochastic  process  F(t)  =  dsflA(s).  Then 
our  first  task  is  to  express  the  stochastic  density  operator 
pst  =  in  the  Kraus-like  operator  representa¬ 

tion  form: 

Pst{t)  =  exp  [~iF(t)az\  p( 0)  exp  [iF(t)az\  (14) 

where  p(0)  =  |^>(0))(V,(0)|  is  the  initial  state  of  the  sys¬ 
tem,  which  is  assumed  to  be  independent  of  the  noise. 

The  desired  Kraus  operators  for  the  spin  are  obtained  by 
taking  a  statistical  mean  over  the  noise  £2^(2)  for  qubit 
A  and  are  given  by 


The  master  equation  for  the  reduced  density  matrix  for 
the  two-spin  system  in  a  non-Markovian  regime  can  be 


Ei  = 


Pa  (l)  0 

0  1 


E, 


qA{t)  0 

0  0 


(15) 
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where  the  time-dependent  Kraus  matrix  elements  are 

QA(t)  =  ^l-p\{t),  and  (16) 

pA(t)  =  exp  [-/(!)],  with  (17) 

f(t)  =  [  G(s)ds 

Jo 

=  + V7t-1)L  (18) 
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where  p*  denotes  the  complex  conjugation  of  p  in  the 
standard  basis  |  +  ,  +),  |  +  ,  —),  \  — ,  +),  |  — ,  — },  and  ay  is  the 
usual  Pauli  matrix  expressed  in  the  same  basis. 

From  the  general  solution  (19-22),  one  can  easily  show 
for  the  initial  state  (23)  that  one  finds 

Q(t)  =  2 max  (\p32{t)\  ~  \J /»n (0)/o44 (0) , 

\pi4{t)  \  -  \/ P22  (0)^33  (0))  .  (26) 


and  similar  expressions  for  pb(0  and  Qb(0-  The  two- 
qubit  case  given  here  can  be  easily  applied  to  N  noninter¬ 
acting  qubits,  an  extension  we  reserve  for  later  attention. 

Since  our  two  spins  are  evolving  independently,  we 
have  the  following  four  Kraus  operators  in  terms  of  the 
tensor  products  of  E\  and  E2: 

*■■  =  ("o  !)® ('»!)■  <19> 

*  =  (»> 

*■•  =  (VS)  «(?!)•  (21) 

a'«  =  ( “  2 )  ®  ( "  2 )  *22) 

IV.  NON-MARKOVIAN  ENTANGLEMENT 
DYNAMICS 

A.  X  matrix  and  concurrence 


B.  Solutions  for  non- Markovian  disentanglement 

The  Ornstein-Uhlenbeck  phase-noise  solutions  for  the 
density  matrix  elements  of  a  general  initial  state  are  given 

by 


Pl'2{t) 

=  Pi2(0)e-M, 

(27) 

Pis{t) 

=  Pi3(0)e-/(t), 

(28) 

P24  (t  ) 

=  *>24(0)e-/(t), 

(29) 

P34{t) 

=  P34(  0)e“/(t) 

(30) 

P23{t) 

=  P23(0)e_2/(t), 

(31) 

Pl4{t) 

=  Pi4(0)e-2^) 

(32) 

Pii{t) 

=  pa(0)  (i=  1,2,  3,  4) 

(33) 

where  f(t)  is  defined  in  (18).  Let  us  note  that  in  the  limit 
7  — >  00,  we  recover  the  standard  Markov  approximation 
where  f(t)  =  Tt/2. 

Although  there  is  no  compact  analytical  expression  for 
the  concurrence  C'(p(t))  with  an  arbitrary  initial  state, 
we  can  readily  show  that  preservation  of  entanglement  is 
restricted  by  the  inequality 


We  now  consider  entanglement  dynamics  of  two  half¬ 
integral  spins  (qubits)  with  an  initial  density  matrix  with 
the  common  X-form  [31]: 


P 


AB 


/  Pll  0  0  P14 

0  P22  P23  0 

0  P32  P33  0 

\P41  0  0  P44, 


(23) 


Such  X  states  occur  in  many  contexts  and  include  pure 
Bell  states  as  well  as  Werner  mixed  states. 

For  two  qubits,  entanglement  can  be  evaluated  un¬ 
ambiguously  via  the  concurrence  function  [33],  which 
may  be  calculated  explicitly  from  the  density  matrix 
pAB .  For  qubits  A  and  B  we  have:  CAB  =  C(pAB  )  = 
max{0,  Q(t)}.  Here  Q(t)  is  defined  as 

Q  =  —  \/^2  —  3  —  \/~^4:  (24) 

where  the  quantities  \  are  the  (generally  time- 
dependent)  eigenvalues  in  decreasing  order  of  the  follow¬ 
ing  (nonlinear  in  p)  matrix: 

C  =  p(aA  ®aB)p*(aA  ®aB),  (25) 


C(p(t))  <  e-2^)C(p(0)).  (34) 

A  sharper  result  occurs  for  the  X  matrix  under  consid¬ 
eration  because  the  diagonal  elements  are  independent  of 
t.  Since  they  all  vanish  as  exp (—/(!))  for  increasing  t,  we 
know  that  Q(t)  must  eventually  become  strictly  negative 
if  diagonal  values  are  initially  non-zero  (e.g.,  any  finite- 
temperature  equilibrium  state).  Negative  Q  mandates 
CAB  =  0,  so  ESD  must  occur.  Next  we  will  consider  key 
limiting  cases. 


C.  Entanglement  decay:  Stationary  limit 

We  consider  now  the  stationary  limit  jt  1.  Then 

/(*)  =  |[*+^(e-7t -!)]-►!*  (35) 

Therefore,  from  (34),  we  get 

C(p(t))  <  e~rtC(p( 0)).  (36) 

Hence,  the  entanglement  decay  rate  is  at  least  as  rapid  as 
T,  and  may  be  much  faster.  Clearly,  the  stationary  limit 
is  identical  to  the  Markov  limit  a(t  —  s)  =  TS(t  —  s). 
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FIG.  2:  The  graph  shows  CAB  vs.  Tt  and  y/T.  The  value 
a  =  1/3  has  been  chosen.  The  reservoir  bandwidth  7  controls 
the  approach  to  the  Markov  limit,  and  we  see  that  for  7  <  F 
the  inevitable  onset  of  ESD  (the  region  where  C(t)  =  0)  can 
be  substantially  delayed. 

D.  Entanglement  decay:  Short-time  limit 

Now  let  us  turn  to  the  opposite  limiting  case:  t<<1. 
In  this  case,  we  can  use  the  following  approximation, 

e-7*  ~  1  -  jt  +  (37) 


where  0  <  a  <  1,  so  the  initial  concurrence  is  £7(0)  = 
2/3[l  -  y/a{\  -  a)]  >  0. 

The  time  dependence  of  the  concurrence  of  this  state 
is  well-known  in  the  Markov  dephasing  limit  [9,  28].  For 
our  present  non-Markovian  case,  which  introduces  the 
environmental  bandwidth  7,  the  time  dependence  of  the 
Q  parameter  satisfies 

=  -<*))■  (42) 

For  all  finite  values  of  7  the  quantity  e-7W  approaches 
zero  exponentially  at  long  times,  and  then  Q(t)  must 
become  negative,  so  ESD  inevitably  occurs,  with  a  finite 
disentanglement  time  £esd  given  by 

e-2/(teD)  =  tt(1_a)i  (43) 

Fig.  3  shows  several  interesting  features  of  entangle¬ 
ment  evolution  under  Ornstein-Uhlenbeck  noise.  Clearly, 
we  see  that  non-Markovian  noises  have  markedly  different 
affects  on  entanglement  evolution  at  short  times,  while 
the  long  time  limit  gives  rise  to  familiar  Markov  behav¬ 
ior.  First,  note  that  ESD  must  occur  for  the  entire  pa¬ 
rameter  range  of  a,  except  for  the  end  points  a  =  0,1, 
given  our  initial  X-state.  However,  the  ESD  times  are  dif¬ 
ferent  for  non-Markovian  short-time  and  stationary  limit 


Therefore, 

m  =  ir  yt2  (38) 

Then  concurrence  decay  is  bounded  by 

C(P(t))  <e~ir^C(p( 0))  (39) 

Hence,  the  effective  disentanglement  time  is  given  by: 


Fdis  — 


(40) 


Clearly,  for  non-Markovian  noises,  the  short-time  limit 
is  more  interesting  since  it  shows  that  the  resultant  en¬ 
tanglement.  behavior  deviates  significantly  from  the  well- 
known  Markov  dynamics.  Obviously,  the  smaller  7  is,  the 
better  approximation  we  have.  For  both  limiting  cases, 
it  is  easy  to  prove  that  a  sufficient  condition  for  ESD  to 
occur  is  pnp22P33p44  ±  0  [34], 


rt 


rt 


V.  EVALUATION  FOR  A  SPECIAL  X  STATE 


In  Fig .  2  we  show  the  evolution  of  the  concurrence  for 
a  specific  X-form  entangled  state: 


pTK  o)  =  - 


fa 

0 

0 

0  \ 

1 

0 

1 

1 

0 

3 

0 

1 

1 

0 

Vo 

0 

0 

1  -  a) 

(41) 


FIG.  3:  The  graphs  show  C'AB  vs.  t,  with  a  =  1/3.  Fig. 
3a  shows  that  with  y/T  =  0.5,  an  initially  non-Markovian 
entanglement  (solid  line)  evolves  in  a  markedly  different  way 
compared  to  its  Markov  limit  (dotted  line).  Clearly,  the  short- 
time  limit  (dashed  line)  gives  a  better  approximation  than  the 
stationary  limit  (dotted  line).  However,  as  shown  in  Fig. 3b 
where  7/r  =  5,  the  difference  is  washed  out  at  later  times 
when  the  short-time  limit  (dashed  line)  ceases  to  be  a  good 
approximation. 
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cases.  From  Fig.  3a,  it  should  be  noted  that  the  dis¬ 
entanglement  time  for  non-Markovian  regimes  could  be 
significantly  longer  than  the  disentanglement  times  in  the 
Markov  limit  if  the  dissipation  is  small.  Flowever,  once 
the  state  becomes  separable  it  will  never  become  entan¬ 
gled  again.  That  is,  entanglement  rebirth  or  revival  does 
not  occur  for  Ornstein-Uhlenbeck  noise  [14]. 

VI.  CONCLUDING  REMARKS 

We  have  presented  here,  as  a  contribution  to  the  sci¬ 
entific  memory  of  Krzystof  Wodkiewicz,  the  first  results 
of  a  new  investigation  that  has  clear  connections  to  his 
long-time  interest  in  quantum  systems  evolving  under  the 
influence  of  stochastic  perturbations.  One  of  the  targets 
of  his  creativity  and  energy,  over  many  years,  was  the 
challenge  presented  by  the  influence  of  non-Markovian 
noise  and  we  have  addressed  that  challenge  with  some 
calculations  focused  on  entanglement. 

We  have  shown  that  entanglement  dynamics  under 
Ornstein-Uhlenbeck  noise  can  be  affected  in  several  dif¬ 
ferent  ways,  depending  on  the  initial  entangled  states  and 
the  noise  correlation  time.  It  can  be  seen  that  the  non- 
Markovian  properties  can  prolong  the  life  of  entangle¬ 


ment.  We  note  that  the  effective  long-time  relaxation 
rate  T  is  ordinarily  associated  with  experimentally  ac¬ 
cessible  relaxation  times  such  as  T\  and  T2.  Fig.  2  high¬ 
lights  the  unusual  domain  7  <  T,  in  which  these  times 
are  shorter  than  the  internal  environmental  relaxation 
time  tc.  Entanglement  survival  is  of  fundamental  inter¬ 
est  at  short  times  in  quantum  information  processing  (see 
[35,  36]).  In  the  case  of  the  short-time  limit  our  results 
capture  the  features  of  quadratic  rather  than  exponen¬ 
tial  decay  at  early  times.  In  this  simple  model,  non- 
Markovian  noises  appear  to  play  a  role  as  a  short-time 
decoherence  buffer,  but  entanglement  measured  by  con¬ 
currence  will  inevitably  conform  to  the  stationary  limit 
at  long  times.  Finally,  it  is  interesting  to  note  that  our 
findings  based  on  the  classical  phase  noise  model  can  be 
extended  into  the  case  of  quantum  phase  noises  where  the 
environment  is  modeled  as  a  set  of  harmonic  oscillators 
at  a  finite  temperature  [37]. 
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